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Abstract 

We investigate the propagation of MHD waves in a homogenous, magnetized plasma in a weakly 
stratified atmosphere, representing hot coronal loops. In most of earlier studies a time-independent equi- 
librium is considered. Here we abandon this restriction and allow the equilibrium to develop as function 
of time. In particular, the background plasma is assumed to be cooling due to thermal conduction. The 
cooling is assumed to be on a time scale greater than the characteristic travel times of the perturbations. 
We investigate the influence of cooling of the background plasma on the properties of magneto-acoustic 
waves. The MHD equations arc reduced to a 1-D system modelling magneto-acoustic modes progressing 
along a dynamically cooling coronal loop. A time dependent dispersion relation which describes the 
propagation of the magneto-acoustic waves is derived by using the WKB theory An analytic solution for 
the time-dependent amplitude of waves is obtained and the method of characteristics is used to find an 
approximate analytical solution. Numerical calculations are applied to the analytically derived solutions 
to obtain further insight into the behavior of the MHD waves in a system with variable, time-dependent 
background. The results show that there is a strong damping of MHD waves that can be linked to the 
widely observed damping of hot coronal loop oscillations. The damping also appears to be independent 
of position along the loop. Studies of MHD wave behaviour in time-dependent background seem to be a 
fundamental and very important next step in developing MHD wave theory applicable to a wide range 
O ■ in solar physics. 
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1. Introduction 

t > , The solar atmosphere is known to be strongly magnetised and dynamic in nature ( |Vaiana et al, 1 973 

|Schrijver et al, 1999 ) . The fine structure of this structured plasma environment consists of a wide range 



O 



in 



C3 



of distinct magnetic features like coronal loops, open flux tubes, prominences, etc. These structures play 
an important role in solar dynamics since they can support various type of MHD waves and oscillations 
which arc thought to contribute to the solution of the long-standing problem of solar coronal heating 
(for some recent reviews see, IKlimchukl l2"00l)l |Erdelyi[ I2DD51 |Taroyan[ I2D051 |Taroyan and Erdelvi[ l2"005|l . 
Observations indicate that waves and oscillations are ubiquitous in the solar atmosphere ( |Wang et al.\ 
gig} [Tomczyk et al] I5U071 iDTPortieu et alA I2UP71 IQkamoto et all I2HTJ71 |Erdelyi and Taroyan[ 121)051 
IVan Doorsselaere et all 120081 Uess et al\ 120091 For recent reviews see, INakariakov and Verwichtel [2005 
|Banerjee ef"ar|l2"0tjTllDe Moortel[IMM|Zaqarashvili and Erdelyi) l2"0"0"9"llMathioudakis et aUl2TjTTj) . Many 
of the above reputs highlight that the waves or oscillations detected are seen to be strongly damped. 
A common suggested mechanisms for damping are resonant absorption ( |Ruderman and Roberts, 2002| 
Goossens et al, 2002| A"schwanden et al, 2003[) and t hermal conduction (Ofman and Aschwandcn, 120021 



Do Moor: i I and Hood 



120031 l2004UDe Moortel eTa/l I2DM1 iMendoza-Briceho et all IMHl |Erdelyi et al\ 



2008P which cause the damping of fast kink and (slow or acoustic) longitudinal waves, respectively. More 



recently, Morton and Erdelyi ((20091 I2010[) argue that radiation of the background plasma may be a 
dominant damping mechanism for coronal oscillations. 

In recent years, coronal waves have been used as a diagonstic tool to infer otherwise unmeasurable or 
hard to measure plasma parameters e.g. magnetic field (fNakaria kov and Ofman, 2001|[Erdelyi and Taroyan, 2008[ )[ 
and scale height ( |Verth et al, 2008"] ). Such a technique is known as magneto-seismology and was suggested 
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by lUchida (1970 ); Zait sev and Stepanov| ([1983")) and lRoberts et a/."l(|1984p in the context of coronal diag- 
nostics. The concept was proposed to be used to any magnetic structure of the Sun by |Erdelyi| (|2006|) 
and labelled as solar magneto-seismology. The models used for deriving the plasma parameters are all 
static in these earlier works, i.e. the background (or equilibrium) plasma is not dynamic, and background 
motion with cooling/heating is more than often neglected. However, there are numerous observations 
which show that the corona, and in particular coronal loops, are highly dynamic in nature, exhibiting 
flows, cooling and heating of the local plasma with a wide range of timescales. For accurate values of 
plasma parameters to be obtained from coronal magneto-seismology this dynamic nature of the plasma 
needs to be incorporated into the models and the importance of a dynamic background must be assessed. 

One feature which should be taken into account when studying the nature of coronal oscillations is 
the temperature evolution of the coronal plasma. There has been a plethora of observations which show 
a number of different temperature evolution scenarios (e.g. Winebarger et al [2003; Naga ta et al. 120031 
|L6pez Fuentes et al.\ I2007P . |Nagata et al.\ (|2003p have observed that there are at least two categories of 
temperature evolution. Hot loops, which are heated to temperatures T > 2.5 MK, are seen by the x-ray 
imagers, e.g., Yohkoh's soft x-ray telescope (SXT) and Hinodc's x-ray telescope (XRT). The hot loops 
are short lived and are seen to undergo relatively fast cooling down to EUV temperatures appearing in 
EUV imagers, e.g., SOHO/EIT. Another category is cool loops which are observed with temperatures 
0.4 — 1.3 MK in EUV images and have relatively long lifetimes. The temperature changes on a slow scale 
when compared to the lifetime of perturbations, or even the loop one's itself and the loops exhibit limited 
dynamic behaviour. The physical process of the plasma cooling depends upon the loop temperature (cool 
or hot). It has been found that radiation is the dominant mechanism for the cooling of the EUV loops 
(T < 2.0 MK) whereas thermal conduction is the cooling method for loops in the region of T > 2.0 
MK. The majority of observed coronal loops that exhibit oscillations are reported to have a temperature 
decrease in an exponential form with cooling times of 500 — 2000 s ( |Aschwanden and Terradas, 2008| 



Ugarte-Urra et al, 20091. 



Propagating plasma disturbances were detected by SOHO/UVCS along coronal plumes (jOfman et all 
HWn [TMfl I2000ap . IDeForest and Gurmanl (TIM5| observed similar compressive disturbances in polar 
plumes by using SOHO/EIT. These observed compressive disturbances are interpreted as propagating, 
slow magneto-acoustic waves, where the suggested damping mechanism for the waves was compressive 
viscosity (Ofman et al, 1999, 2000b). This problem was further investigated bv iMcndoza-Bricch o et al\ 
(|2004p who considered a gravitational stratification in addition to dissipative processes of thermal con- 
duction, viscosity and radiation. They found that oscillations are even more efficiently damped in stratified 
plasmas. Similar intensity disturbances were observed in coronal loops by TRACE ([Nightingal e et al, 1999| [ 
Schrijver et al, 1999||De Moortel et al, 2000||McEwan and De Moortel, 2006[ ) and EIT/SOHO perghmans and Clet 
Nakariakov et al] (|2000p have also identified these disturbances as slow magneto-acoustic waves and 
they found that the wave evolution is affected by dissipation and gravitational stratification. Moreover, 
|Wang et aj]$M5ty and |Taroyan et aZ.|([2"0tJ7) have observed (T > 6 MK) hot loop oscillations by SUMER 
and SXT. They suggested that these oscillations are standing longitudinal (acoustic) waves and are 
triggered by a footpoint microflare. It was reported that after an initial temperature increase the loop 
is observed to be cooling. Recently, |Erdelyi and Taroyan (|2008|) and Wang et al. (|2009p have detected 5 
minutes quasi-periodic oscillations in transition region and additional five coronal lines by Hinode/EIS 
at the footpoint of a coronal loop. These oscillations were interpreted as slow magneto-acoustic waves 
propagating upward from the transition region into the corona. Further, it was found that the amplitude 
of oscillations decreases with increasing height. It is suggested that the source of these oscillations 
are the leakage of the photospheric p-modes through the chromosphere and transition region into the 
corona, i.e. a similar mechanism that is put forward to explain spicule formation ( De Pontieu et al, 2004[ 
|De Pontieu and E rdelyi, 2 006| ) and transition region moss oscillations ( |De Pontieu et al, 2003 ). 

A theoretical study of how the slow mode under solar coronal conditions can be extracted from the 
MHD equations has been done by, e.g. [Roberts (2006]). The magnetic field was assumed straight in the 
vertical direction in a gravitationally stratified medium. This theory has been applied to the slow waves in 
coronal loops observed by SUMER and TRACE. ILuna-Cardozo et al.\ (|201ip generalised these studies to 
concurrently occurring magnetic and density stratification, i.e. taking into account the strongly expanding 
nature of loops that support MHD waves. 

The hot loops are thought to be heated by impulsive heating events, e.g. flares, micro-flares, nano- 
flares. The origin of these energisation mechanisms could be either waves or reconnection (jAntolin et al\ 
I2008al I2008bl I2009p . There have been a number of simulations investigating the evolution of coronal 
loops after heating events (see, e.g. Uakimiec et aTX 119921 [Cargill[ H994UAschwanden and Tsiklauril 120051 
Taroyan and Erdelyi) [20091 [20T01 [Taroyan et al.\MT\\<i . Impuls lve heating events only last for a relatively 



short duration and if the plasma reaches a hot enough temperature by the end of the heating, thermal 
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conduction dominates the cooling of the plasma. From calculations by e.g. |Cargih] (|1994p . the decrease 
in temperature due to thermal conduction takes an almost exponential form. 

The effect of thermal conduction on MHD waves has been extensively studied. However, all the pervious 
models have assumed the background plasma is static, i.e. no change in the equilibrium plasma quantities 
over time. Ilbanez and E scalona (|1993[) studied the propagation of thermal and magnetosonic waves in 
optically thin plasma. They found that thermal waves are always damped whereas the magnetosonic 
waves are damped in the range of temperature 10 4 K < T < 10 s K due to thermal conduction mecha- 
nism. Moreover, |Ofman and Wang] (|2002p and IMendoza-Briceno et al.\ (|2004p have used a ID nonlinear 
dissipative MHD model to study the Doppler shift oscillations of hot coronal loops observed by SUMER. 
The oscillations were interpreted as slow magnetosonic waves. Compressive viscosity was found to be less 
significant on damping loop oscillations compared to the thermal conduction mechanism. The propagation 
of MHD slow waves in a static ID isothermal medium has also been studied by IDe Moortel and Hoodl 
(|2003p where thermal conduction and compressive viscosity are considered as damping mechanisms for 
perturbations. They found that 'enhanced' thermal conduction is the dominant damping mechanism of 
slow waves in coronal loops. Further to this, IDe Moortel and Hoodl (|2004[) investigate the behaviour of 
slow MHD waves in a stratified and diverging static atmosphere. The introduction of stratification reduced 
the effectiveness of thermal conduction on the damping of the slow wave. This result is in disagreement 
with those of Mendoz a-Briceho et al.\ (|2004[) . who found an opposite effect for slow standing waves. 
Further, there is reported a decay of oscillation amplitude by optically thin radiation and area divergence 
although this is only found to change the rate of damping by a few percent. Overall, it is found in 
many earlier studies that thermal conduction is an essential cause of the damping of loop oscillations in 
static equilibria when compared to the other mechanisms under the coronal conditions (i.e. T = 10 6 K). 
However, there is an important point to highlight: although thermal conduction seems to be an efficient 
dissipative mechanism among the various damping mechanisms, considered, it can not efficiently model 
the observed decay of amplitudes, when coronal parameters are substituted into the classical formulae of 
thermal conduction obtained by |Braginskii| (|f 965|) and |Chapman and Cowling] (|1958"| . etc. This means 
that modelling efforts need to consider alternatives. We provide here one. 

The damping of oscillations due to a cooling background coronal plasma is a novel and physically 
natural idea. Morton and Erdelyi| (|2009l) found that the cooling of the background plasma could lead 



to the damping of even fast kink (i.e. by largely incompressible) oscillations. In addition to this, they 
noticed that the ratio of fundamental mode to the first harmonic decreases as the loop cools. Furthermore, 
|Morton and Erdelyi| (|20I0p found that the damping due to the cooling of the plasma could account for 
the observed damping in a number of examples of observed transverse oscillations. At the time of writing, 
this number was almost all the known examples available in the literature (|As chwandc n et ai[ 119991 120021 
INakariakov et ai[ I1999[) . In a recent theoretical work by IMorton et al. \ (|2010[) , the effect of a radiative 
cooling background plasma state on the propagation of magneto-acoustic waves in a uniformly magnetized 
plasma was investigated. Although the approximation of unbounded uniformity of the plasma may seem 
to be a severe simplification, this step was necessary in order to give insight into the underlying physics. 
The radiation mechanism was assumed to be the mechanism for the cooling of the plasma and the plasma 
cooled exponentially in time, which is capturing well the key features of the observational data. As a result, 
they found that the slow and fast modes are damped due to cooling. The radiative cooling was shown to 
damp the slow mode by up to 60% within characteristic lifetimes. This effect is much stronger than the 
contribution predicted by IDe Moortel and Hoodl (|2004p and there was no need of additional, less realistic 
assumption about the coefficient of thermal conduction. 

In this paper, we study further the propagation of longitudinal MHD waves in a uniform magnetized 
plasma under the influence of a cooling background state. The cooling of the plasma is assumed, here, 
to be dominated by thermal conduction, so is applicable to oscillations in hot loops. Solutions to the 
background plasma equations show the cooling profile of the plasma is exponential in time in agreement 
with observational reports. The dispersion relation which describes the slow modes and their properties 
is derived. An analytic expression for the time-dependent amplitude of the waves is also derived. The 
results show that thermal conduction has a dominant effect on the slow modes where their amplitude 
is strongly damped in plasmas with a dynamic background state. The rate of damping is quantified in 
the model and is found to be dependent upon the amount of stratification in the plasma and the initial 
temperature of the background. The damping is also shown to be weakly dependent upon position along 
the slowly varying equilibrium. 

2. Governing equations 

Consider a homogeneous magnetised plasma where the background temperature is changing as a function 
of time due to thermal conduction and density is constant in time. The magnetic field is uniform and 



SOLA: latex-attempt.tex; 12 November 2010; 1:29; p. 3 



R. Erdelyi, K.S. Al-Ghafri and R.J. Morton 



in the z direction, i.e. Bo = BqZ. The governing MHD equations for the background plasma take the 
following form 

^+V.(pv)=0, (1) 

<9v 1 

+p(v.V)v = -Vp+ — (V x B) x B + pg, (2) 
at fiQ 



R pi 
A (7 - 1) 



d T . T 

- + (v.V) 



9B 

~dt 



k\7 2 T, (3) 

p> -j 

— =Vx(vxB), (4) 

P = ypT, (5) 
P 



where v is the flow velocity, B is the magnetic field, g is the gravity, po is the magnetic permeability 
of free space, 7 is the ratio of specific heats, R is the gas constant, p, is the mean molecular weight, T 
is the temperature, kV 2 T is thermal conduction term, p and p are the plasma density and its pressure, 
respectively. 

The medium is assumed to be cooling due to thermal conduction with no temporal change in density 
and it is also assumed that there is no background flow (i.e. vo = 0), so Eqs. (JTJ — (J5} reduce to 

vo=0, ^=0, (6) 

Vpo = Pog, (7) 

Po = po(z), (8) 
R 

Pa = —PoTq, (9) 
P 

m - kW 2 T , (10) 



£(7-1) dt 



where the (0) index denotes background quantity. 
Since Eq. can be written as 

—dp = ^dz, (11) 
Pa H 

where z = Iz (here / is some characteristic length scale) and H = po/pog, and we are interested in a small 
characteristic scale length which is much smaller than the scale height, H, then we arrive at 

Vp w 0, (12) 

which is representative to a weakly stratified atmosphere, i.e. l/H <C 1. The solution to Eq. (TTOl) by 
separation of variables, and using Eq. ((9j> gives the temperature as 

T (z,t) =T ie xp [ - (7 Z^M. ) [_ C2Z 2 + C3Z + C4 ] ] (13) 



Rpi 



and the density as 



p (z) = 2 , (14) 

-C 2 2^ + C 3 Z + C 4 



which is physically valid when €2,03 -C 1 and C3 > C2, where A is the separation constant, C2 = A/2, C3 
and C4 are constants, Ti is the initial temperature at z = and i = 0, and pi is the density at z = 0. 

Note that the background temperature decreases exponentially with time. It has been shown for radia- 
tive cooling loops that an exponential profile provides a good fit for the observed cooling ((Xschwandcn and Tcrra dasl [ 
120051 |Morton and Erdelyi) 120051 12010)1 . To the best of our knowledge, there are no such observations for 
hot loops cooling due to thermal conduction, that is why dpo/dt = here. 

The form of the density profile in Eq. ([T4")) is consistent with a weakly stratified atmosphere. For a 
gravitational stratified atmosphere the density profile would be formed to be 

po = piexp ■ (15) 
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The hydrostatic scale height H is given by H = 47T 



Mm 



MK 



hence for hot loops H is in the range 94 — 282 



Mm. In the case of weak stratification, i.e. l/H <C 1, the density can be approximated as 



Po 



z 1 , 
l-r= + r — 



H 2 \H 



(16) 



Eq. (fT4]) can be written 



Po = 



1 



pi(l -y)~pi[l- (~c 2 z 2 + c 3 z + c 4 - 1)] , 



upon the assumption y is small. Comparing to Eq. (|16|) . then the values for the constants are C4 = 1, 

C3 = jj and C2 = -| = \ (jj) 2 and hence ?/ is small. 

Perturbing the background equations, the linearized MHD equations can be found by writing all 
variables in the form 

f(z,t) = f (z,t) + f 1 (z,t), 

and neglecting all terms containing squared perturbed variables. Here the subscript (0) represents the 
equilibrium quantities and the subscript (1) indicates the perturbed quantities. 

Since thermal conduction is known to have a strong effect on slow modes, we will concentrate our 
analysis on the properties of slow modes. It has been shown by a number of authors, e.g. IRobertsI (|2006p 
and lLuna-Cardozo et aZ.I(|2011[) . that the slow modes can be isolated by assuming v\ x = v\ y = 0. Therefore 
the linear dissipativc MHD equations for the slow modes in the presence of thermal conduction reduce 
to a 1-D system given by 



dpi 
dt 



Po- 



dvi 



dpo n 
oz 



Po- 



dvi 
dt 



R 
P 



Pi dT po dTi 
7-1 dt 7 - 1 dt 
dBi 



PqTq 



dvi 
dz 



To dp 

7 — 1 oz 



dpi 

dz ' 

nV 2 Ti, 



= V x ( Vl x B ) = B ^ - B (^i), 
at oz oz 

R 

Pi = — [poTi + piT ] , 
P 

V-Bi =0. 



(17) 
(18) 
(19) 

(20) 

(21) 
(22) 



Here vi = vi z . Dimensionless variables can now be introduced to simplify the equations, where the 
following dimensionless quantities are suggested: 



Pi 



Pi 



Po 

IVi 
Pi 



Po 



Pi 



Pi 



, Ti 



Ti 



, T 



To - 



Vl 



(23) 



t = — , z 

T 



I 



H 



62 = lC2, C 3 = lc 3 . 



Here c S i is the initial sound speed, / is the wavelength of the oscillations and r is the sound travel time 
of a wavelength. Thus the equation of energy in terms of dimensionless variables removing tilde, will be 



dTo dTi dp 
P 1 ~dt + f ) °~dt + ^~ 1 )f )oToW - Vl ~ T ° Vl ~dz~ 



{ d 2 Ti 
dz 2 



where d 



(7 - l)npiTi 



(24) 



Using coronal values, (see e.g. IDe Moortel and Hood! I2003[) . we find d is a small quantity, where the 
standard coronal values of all variables 



' To = 


= 1-6 MK, 




Po = 


= 1.67 x 10~ 12 kg 


m" 3 , 


K = 


10 _H T 5/2 w m - 


- 1 deg" 1 


« P = 


0.6, 




R = 


8.3 x 10 3 m 2 s" 2 


deg -1 , 


7 = 


5/3, 




, T = 


300 s, 





(25) 
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give a value of d = 0.04 for T = 1 MK and d = 0.61 for T = 6 MK. 
From the linearized ideal gas law equation 



Px = PoTi + piT , 



(26) 



and the continuity equation 



dpi 

at 



Po 



dvi 
dz 



Vl 



dp 

dz 



(27) 



the energy equation takes the following form 



dpi . „ vvx ,< 
-xr +7Po-7o— = d- . 
ot Oz oz z 

Differentiating Eq. (|28l) with respect to z and using Eq. (|T8|l . we arrive at 



dv x ,d 2 Tx 



d 2 vx 
~dt 2 ~ 



d 2 vi 
~dz T 



d d 3 Tx 
IPo dz 3 ' 



(28) 



(29) 



To solve the equation of energy, we start from the continuity equation and using the gas law equation to 
find an equation in terms of velocity variable. 
Substituting Eq. into Eq. ([27|). we obtain 



1 dpi _ px dT 
T dt T 2 dt 



Po dT x i p Tx dT 
T dt 



T 2 dt 



Po 



dvx 
dz~ 



The dimensionless background temperature is 



Mz, t) = 



exp(-Xdt) 
Po(z) 



Vx 



dpo 

dz 



0. 



(30) 



so 



dTo 
dt 



= STq, where 5 = —Ad. 



Then, after substituting ([3"Tj) and multiplying by To, Eq. (|3T)|) will be 



dpx x u±i uvx upo 

— - dpx - Po^tt + OpoTx + poT -^- + T vi—- 

ot ot oz oz 



dTx 



dvx 



dpo 



0. 



Substituting this equation into Eq 

d 2 vx 



, we obtain the following equation 



dvx d dp d 3 Tx 

(7 1)poTo dW- z w ¥ - dl {ToVl) H = d dldz~ 2 



Po- 



d 2 Tx 
dt 2 



5p °-dt> 



(31) 



(32) 



(33) 



which represents another relation between the perturbed temperature and perturbed velocity in addition 
to Eq. ([29]) . Now, we aim to reach the governing equation for the velocity perturbation. Differentiating 
Eq. (1321) once with respect to t and three times with respect to z, and using Eqs. (j28|) and ([29)) with some 
cumbersome algebraic operations, we arrive at 



7' 



Po d 4 vx 

d at 4 



+ 7 

j6p T d 2 



Po& 



35 ( dp { 



dpo 
Po V dz 



2 d 2 po 1 dp d 



i:) 2 



Po dz 2 po dz dz dz 2 



d 3 



i ; i 



P 2 o\d 



25d 2 p 6 dpo d | pqTq d 2 
po dz 2 po dz dz d dz 2 



dz 2 



dz dz 3 dz A 



dvx 
~dt 



dT Q d 3 d 4 
Oz oz 6 Oz* 



dt 3 
d 2 vx 

%ix = 0. 



(34) 



To simplify Eq. (|3"4"]) we neglect all terms of order dS and higher degree in d or 6 as d and S are small 
parameters. Then, after multiplying equation (|34[) by d, we obtain the governing equation 



7Po- 



d 4 vx 
W 



IPoTo 



d 2 d 2 Vx 

dz 2 dt 2 



3d f dpo 
pod + -o 

Po 



dz 



2d d 2 p 
Po dz 2 



d dpo d ^d 2 
po dz dz dz 2 



d :i 



vx 



dt 3 



a 4 



d 2 dT d 3 
jSp T ^-2d—— 3 -dT — i 



dvx 
~dt 



(35) 
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The background quantities are slowly varying when compared to the perturbed quantities, so eliminating 
the small terms gives 



d_ 

at 



d 2 v\ _ d 2 V! 
dt 2 dz 2 



= 5 



d 2 v\ 



d 2 vi 
Ih 2 



d__^_ 

Po dz 2 



d 2 v 1 T d 2 v x 



dt 2 



7 dz 2 



(36) 



In the case of no thermal conduction, Eq. (|36[) reduces to the wave equation that represents (longitudinal) 
acoustic mode propagating in a flux tube and has the form 



d 2 vi 

~dt 2 ~ 



t d2vi - n 
dz 1 



(37) 



while in the case of unstratified atmosphere, i.e. no change in the background plasma quantities, Eq. (1361) 
leads to the model governing equation found by IDe Moortel and Hoodl (|2003p as follows 



d_ 

at 



r d 2 v\ 



dt 2 



T 



d 2 vi 



dz 2 



dz 2 



d 2 v\ 1 d 2 v\ 



dt 2 7 dz 2 



(38) 



Eq. (|36[) can now be checked. Assuming the background is constant in time and space, i.e. static and 
unstratified, one can Fourier-analyse all perturbations, exp(i(wi — kz)). Hence, Eq. ([55)1 reduces to the 
dispersion relation 



idui 2 k 2 



L0k 2 +i-k A 



7 



0, 



(39) 



where uj is the frequency and k is the wavenumber. This dispersion relation has been obtained first by 
IFieldl (|1965|) . This provides confidence that Eq. (l36l) is consistent with earlier studies (see, e.g. IFieldl 119651 
and IDe Moortel and Hoodl 1MB)) . Note that, in the limit d — > 0, i.e. no thermal conduction, Eq. (|36D 
reduces to Eq. (|3"T)) i.e. uj 2 = k 2 =^> uj = ±fc which is, rightly, the undamped (longitudinal) acoustic 
mode. 



3. Analytical solutions 

We now seek to find an analytic solution to the governing Eq. (|36p . Unlike previous models which 
include thermal conduction, there is now a time dependence due to the temporally changing background 
temperature. This means we can not Fourier analyse in time. Instead, since the RHS of Eq. has 
derivatives multiplied by a small factor, d, the WKB approximation (see, e.g. |Bender and Orszag| I1978P 
can be applied to find an approximate solution. The WKB approximation has a greater accuracy for the 
smaller value of d, however, it can still provide an accurate solution even for d « 1. Since the variables in 
Eq. (|36[) depend on time t and space z, we introduce 'slow' variable t\ = dt and 'local' variable Q — dz to 
solve this equation. The slow timcscale physically means that the conductive cooling timescalc is longer 
than the period of the oscillations. The 'local' lengthscale means any significant changes in quantities in 
the z direction occur on length scales longer than the wavelength of the oscillation. The velocity in terms 
of the new scaled variables and by the WKB approximation has the form 

wi(C, *i) = Q(C, exp U9(C, h)j , (40) 

where Q((, ti) and 0(£, ii) arc functions to be calculated. 

Substituting Eq. (|40)l into Eq. ((36)) and taking the largest order terms in d, which is d -4 , we obtain 

(j|) 2 - c 2 s (z, t) {^j = 0, c s (z, t) = Vn. (41) 

If we assume d&/dt = co(z, t) and dO/d( = k(z, t), where u>(z, t) is frequency and k(z, t) is wavenumber, 
then Eq. (|4ip is in fact a temporally and spatially dependent dispersion relation for the longitudinal 
(acoustic) mode. 

The next largest order terms in d (of order d~ 3 ) give the equation for the amplitude, 

dQdQ 7, BQ dc s BQ , „{dQ\ 3 n 2 dQ dQ n , ln . 

- 2lPoCs W 1 ^ + 2 XpoCsQ W~ 1P ° Cs ^ Q ^~ h ~ )CsQ {~d() +2Ws 9C9C ' ( } 
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which, after some algebra, has the form 



-1 dQ dQ 



4^ 



i dc s (7-1) fee 



2c s d( 27p c s V 



O = o. 



(43) 



Eqs. (f¥T|) and (f4"3")l will be solved by using the method of characteristics so we need to derive boundary 
conditions at z = 0. 

To achieve this we study a thin layer around z = where we may assume the spatial gradients of both 
To and po are very small, so they can be considered constant in space in this region. This enables the use 
of Fourier-analysis such that, there is ^exp(ikz), for perturbed variables. Eq. (|5o) reduces to 



IPo- 



"f[p$\d + dk' 



Vl 



dt 3 



IPoTok 



,d 2 vi 
"di? 



- hXdp Q T Q k 2 - dTofc 4 ]^ = 
dt 



(44) 



Next, we apply the WKB method to this equation by assuming that the perturbation in terms of the 
variable t\, where t\ = dt, has the form 



)i(ti)cx P ( ^ei(*i; 



Substituting into Eq. (|4"4")) . the highest order equation in d gives 

dQi 



dti 



Eq. (|4"6")l has the solution 



The next order equation in d is 



which has the solution 



= c s (ti)k. 



2k, 



= y [l-c fl (ti)]. 



dQi 



A _ ( 7 - 1) k 2 

i 7 y 



Qi(£i) = cxp 



A (7-l)fc 5 



= 0, 



ti 



7 



(45) 



(46) 



(47) 



(48) 



(49) 



where k is the wavenumber at z = 0. 

We now have the necessary approximate boundary conditions at z = to proceed with the solutions 
of Eqs. (|4"Tj) and (j43]) using the method of characteristics. 

Then, the solutions of Eqs. (|4"Tj) and (|4"3")l can be found by the method of characteristics and we 
introduce the variables r and s to parameterise these equations. The characteristic equations of Eq. (1411) 
are given by 

dh -l 9c = ae 

s c s os os 

with boundary conditions for the characteristics curves 

h(r,0) = r, C(r,0) = 0, 0(r, 0) = 0i(r), 

where 0i(r) is given by (|T7|) and (r, 0, 0i) is a point in the curve (ti(r, s), C(r, s), 0(r, s)) which lies 
in the surface {(ii, (, ®((,ti))} that represents the graph of the function 0(C, ii). The solutions of the 
characteristic equations are 



s(C,ii) = C, 



r(C,ti) = -^ln 



.—Xti. Xd 
exp( ) — - — = arctan 



e(c,*i) = 



-2k 



ex P(— 



arctan 



c 3 7fxo+2c 2 (c-if) 
k 2vW^xo-<v^(c-if) 

c 3 /F(C) + 2 C2 (C-f^) \ 
2^r 2 ^F%)~c^ 2 (C~^)) 



- 1 



(50) 
(51) 

(52) 
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Similarly, the characteristic equations of Eq. (|4"3")l are 
-1 



~ds~ 



with boundary conditions 



ds 



i. 



dQ 
ds 



A 

4c., 



i dc s (7-1) /ae 



2c s dC %lPoc s \ dC, 



Q- 



h(r,0)=r, C(r,0) = 0, Q(r, 0) = Q x (r), 

where Qi(r) is given by (|4T))) and (r, 0,Qi) is a point in the curve (ti(r,s),(( r ,s),Q(r,s)) which lies 
in the surface {(£1, Ci Q(C)£i))} that represents the graph of the function Q(C, £i). The solutions of the 
characteristic equations are Eqs. ([50|) . (|5Tj) and 



Q(C,£i) = exp 



7 — 1 k 2 



d , Xt\ . 
exp( ) arctan 

a/C2 2 



c 3 7f(C) + 2c 2 (C^if) 

2VWi^-C3^(C-if) 



111 



*X0 



where 



^(0 



-c 2 C + c 3 dC + <r 



MC) 



(53) 



(54) 



Taking Taylor expansion of the amplitude and ignoring the small terms under the assumption of slowly 
varying background coronal plasma and weakly stratified atmosphere, we arrive at 



Q(C,£i) = i 



7 - 1 k 2 



(55) 



Eq. (|4T)]) . combined with Eqs. ([5~Tj) — (|54")l . provides a full solution to Eq. All other perturbed quantities 
can now be calculated using Eqs. (|17l) — (|22|) . The equations may not reveal much information yet as they 
are very complicated in nature. We examine a limiting case for the amplitude to reveal some properties. 
Assume that A = (unstratificd atmosphere), then the system of equations (|5Tj) . (|53|) and ((54)) reduces 
to 

( 7 -l)fc 2 ' 



QiCh) = exp 



7 



(56) 



This limit is in good agreement with its counterpart in lDe Moortel and Hoodl (|2003 |1. The amplitude of 
the wave is damped and the damping is dependent upon the measure of d and k. 



4. Numerical calculations 

It has been shown in IMorton et al. \ (|2010p that the WKB estimates provide good approximations to the 
frequency and amplitude variations in time and space when the plasma is cooling due to radiation. We 
see no reasons why this should be any different here, where thermal dissipation in the format of thermal 
conduction is considered. Numerical calculations are now used to evaluate how the analytic solutions of 
the MHD waves behave in a system with variable background. 

The amplitude of slow wave is computed using characteristic coronal values. In Figure [1] we plot the 
amplitude of the longitudinal (acoustic) wave for different values of A (which is the separation constant 
and is related to the squared reciprocal of scale height) and the value of thermal ratio, d, at z = 0. It is 
found that the amplitude decays rapidly in time for all values of A and d, over typical observed timescales 
of oscillations. In particular, Figure [Tal depicts that the oscillation amplitude declines dramatically under 
the coronal values of temperature 10 6 K. Figure [lb] illustrates a faster decay of the wave amplitude in 
a region of temperature 3 x 10 6 K, i.e. for even hotter loops. The decay of the amplitude is enhanced 
even further for very hot (SXT/XRT) loops, i.e. T = 6 x 10 6 K as evidenced by Figure [Tel In each plot 
shown in Figure [TJ it is also clearly seen that increasing the stratification of the plasma, i.e. increasing 
A, decreases the rate of the damping. This is consistent with lDe Moortel and Hoodl (|2004[) . who showed 
that adding stratification decreased the rate of damping of the slow mode due to thermal conduction. 
The new and important result here is that the rate of damping is much stronger than that predicted by 
IDe Moortel and Hoodl (|2004|) . Further, there is no need here to artificially amplify the heat conduction as 



SOLA: latex-attempt.tex; 12 November 2010; 1:29; p. 9 



R. Erdelyi, K.S. Al-Ghafri and R.J. Morton 



opposed to lDe Moortel and Hoodl (|2004|) . A note of caution is, however, needed when interpreting these 
results. The temperature of the plasma drops rapidly due to the cooling. As T — > 2 MK radiation will 
start to become the dominant damping mechanism. So as T reaches 2 MK the damping of the longitudinal 
mode may not look as shown here. A full numerical simulation incorporating the variable background and 
the variable value of d will need to undertaken that is well beyond the scope of the present study and may 
require a full numerical approach even for an initial insight because of the far to complex mathematical 
modelling (e.g. generating background flows, etc.). 

In Figure^] we show how the magnitude of the thermal conduction coefficient, k, effects the rate of 
damping of the longitudinal (acoustic) mode. Altering the value of k by an order of magnitude leads to 
large changes in the rate of damping, i.e the amplitude of longitudinal (acoustic) wave decreases rapidly 
with an increasing value of n. 

Next, Figure [3] shows how the decay of the amplitude varies for different values of z = (0.0,0.1,0.5). 
Two values of thermal ratio d = (0.04, 0.22) are taken to illustrate the influence of thermal conduction on 
cool and hot corona, T=(l MK, 3 MK), and the characteristic value of A = 0.1. This plot shows that the 
strength of the damping due to thermal conduction will be for all intents and purposes, the same along 
a loop. 



Oscillations Amplitude Oscillations Amplitude Oscillations Amplitude 




300 600 900 1200 1500 300 600 900 1200 1500 300 600 900 1200 1500 



t (s) t (s) t (s) 

Figure 1. The graphs show the amplitude of oscillations with different values of A (0.1, 0.5, 0.7) characterising stratification 
and specific value of d, i.e. the value of thermal ratio, at z = 0. (a) d = 0.04 (T = 1 MK), (b) d = 0.22 (T = 3 MK), (c) 
d = 0.61 (T = 6 MK). 



5. Discussion and Conclusion 

The influence of variable background on the longitudinal (field-aligned) MHD wave propagating in a 
homogeneous, magnetised plasma in a weakly stratified atmosphere has been investigated. The back- 
ground plasma is assumed to be cooling as the magneto-acoustic wave propagates. Thermal conduction 
is the dominant mechanism causing the plasma cooling and changing pressure as a function of time. The 
magnetic field was assumed to be constant and directed along the z direction. A governing equation is 
derived under the assumption that the background plasma is cooling on a time scale greater than or 
comparable to the characteristic period of the perturbations. A time- and space-dependent dispersion 
relation that describes the propagation of the longitudinal magneto-acoustic mode is obtained and an 
analytic equation for the time-dependent amplitude is also derived. The effect of cooling background due 
to thermal conduction on the amplitude of hot loop oscillations was then studied. 

The governing equation in terms of time and space is solved by using the WKB theory. Leading and 
first order equations for the dispersion relation and wave amplitude, respectively, are obtained and solved 
analytically. An approximate solution which represents the properties of the field-aligned acoustic wave 
is found with the aid of method of characteristics. Numerical evaluations illustrate the behaviour of the 
analytic solutions. A comparison of wave behaviour for a range of initial temperatures is studied. 

Enhanced by the cooling background plasma, the amplitude of longitudinal (acoustic) waves was found 
to rapidly decay due to the influence of thermal conduction. The rate of damping of the oscillations 
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Oscillations Amplitude 




t(s) 

Figure 2. The graph shows the amplitude of oscillations with different values of the thermal conduction coefficient, 
K = (lCT 10 , fO^SlCr 12 ) at z = and A = 0.1 where T = 3 MK. 

was found to depend on the initial temperature of the plasma and the amount of stratification. It was 
previously shown that thermal conduction leads to the damping of slow mode oscillations. For example, 
the amplitude of slow modes in IDe Moortel a nd Hood ( 20Q3]) was detected to undergo damping whereas 
in this study the amplitude was found to experience a much stronger damping. It should be noted that 
IMorton et al.\ (|2010p reported that the damping of the slow mode was also heavy where the amplitude 
of wave suffers a rapid decrease due to the cooling of the background plasma by another dissipation, 
namely by radiation mechanism. Consequently, we conclude that the magneto-acoustic oscillations of hot 
corona can be exposed to a strong damping because of cooling by thermal conduction while radiation is 
a dominant method of damping cool coronal oscillations. 

Overall, the presented results into the study of the damping of coronal oscillation imply that both 
radiation and thermal conduction mechanisms should be taken into account even for background plasma 
as a result of their efficient effect on dominating the properties and lifetimes of MHD oscillations in 
cool and hot corona, respectively. The main point is that these dissipative processes introduce a dynamic 
plasma background. Such background in often reported in observations. The MHD wave theory of dynamic 
plasma is in its infant stage and we argue that for an adequate modelling of the solar coronal processes 
such studies as consideration of dynamic plasma background in the modelling is essential. The temporal 
evolution of the background plasma has an extremely important effect on the properties of waves and 
needs to be incorporated into future models. This is especially necessary if deriving plasma parameters 
from magneto-scismological methods, e.g. in applications to the solar corona. 
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Figure 3. The graph shows the amplitude of oscillations as function of time at variable positions along the coronal magnetic 
field lines, e.g. z = (0.0, 0.1, 0.5) and A = 0.1, and with different values of d = (0.04, 0.22). 
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